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Abstract. Critical behavior developed near a quantum phase transition, interesting 
in its own right, offers exciting opportunities to explore the universality of strongly- 
correlated systems near the ground state. Cold atoms in optical lattices, in particular, 
represent a paradigmatic system, for which the quantum phase transition between 
the superfluid and Mott insulator states can be externally induced by tuning the 
microscopic parameters. In this paper, we describe our approach to study quantum 
criticality of cesium atoms in a two-dimensional lattice based on in situ density 
measurements. Our research agenda involves testing critical scaling of thermodynamic 
observables and extracting transport properties in the quantum critical regime. 
We present and discuss experimental progress on both fronts. In particular, the 
thermodynamic measurement suggests that the equation of state near the critical point 
follows the predicted scaling law at low temperatures. 



PACS numbers: 64.70.Tg, 67.85.Hj, 64.70.qj 
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1. Introduction 

A quantum phase transition occurs in a many-body system at absolute zero temperature 
when a new order smoothly emerges in its ground state [U |2] . Even though a quantum 
phase transition happens at zero temperature, the transition influences the flnite- 
temperature properties of the many-body system, which leads to possible experimental 
observations. Near a quantum phase transition, the many-body system shows universal 
scaling behaviors, termed as quantum criticality, with characteristic scaling exponents 
determined by basic properties of the system, such as symmetry and dimensionality 
[1]. Quantum criticality is of wide research interest because both spatial and temporal 
correlations in such systems are expected to follow intriguing universal relations. 
Understanding quantum criticality represents a major challenge to our knowledge of 
many-body physics. 

In contrast to classical phase transitions driven by the competition between energy 
and entropy of a system, quantum phase transitions are driven by the competition 
between different quantum energy scales [11 [3l H] . In a system of bosonic atoms confined 
in optical lattices, the superfluid-to-Mott insulator transition, described by the Bose- 
Hubbard model [3], is ideal for studying quantum criticality. In this system, the two 
competing energy scales are the tunneling energy t and the on-site interaction U . At 
zero temperature, two phases are predicted: the superfluid and Mott insulator phases, as 
shown in Fig. [H In the tunneling-dominated regime {g = t/U ^ 1), the system forms 
a superfluid at low temperatures; in the interaction-dominated regime {g ^ 1), the 
system assumes a Mott insulator phase in the ground state. When g equals the critical 
value gc, the ground state of the system is neither a superfluid nor a Mott insulator. At 
finite temperatures, the critical point expands into a V-shaped regime where universal 
scaling behaviors are expected [H O Ej , see Fig{Tl 

The flexible control of cold atoms and the lattice potential provide essential tools 
to study quantum phase transitions and quantum criticality. In optical lattices, the two 
competing energy scales, tunneling and on-site interaction energies, can be controlled 
precisely by tuning parameters of optical lattices or the magnetic fleld near a Feshbach 
resonance [7] . Using evaporative cooling [HI |9] and other techniques [10], [11] , one can 
reach a temperature low enough that the thermal fluctuations are largely suppressed. 
Furthermore, high resolution imaging techniques have been recently developed to 
measure the density proflles and fluctuations of the trapped atomic gas P [HI [12]. 

In addition to the optical lattice, atoms are typically conflned by an external 
harmonic potential, which adds new perspectives for experimental observations. Because 
of this potential, the system has a higher density as one moves toward the trap center 
[3], [12]. This density proflle can reveal the phase diagram of a homogeneous system 
with different chemical potentials at a flxed coupling constant g, as shown in Fig. 
M^a). A system with a large and negative chemical potential is in a vacuum state. As 
one increases the chemical potential, the system evolves into the superfluid state. For 
small enough g and larger chemical potential, the system can enter the Mott insulator 
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state with unit occupation number. Typically, the harmonic potential is slowly varying 
compared to the optical lattice potential, and the density distribution is expected to 
be slowly varying over the length scale of the lattice constant. With a weak harmonic 
confinement, the local density approximation applies, and every point in the system can 
be viewed as a homogeneous sub-system with a local chemical potential. Therefore, a 
measured density profile probes a line in the phase diagram along the chemical potential 
direction. Near the phase boundary, the line crosses the quantum critical regime of our 
research interest, as shown in Fig. |2l 

The main goal of this paper is to describe two possible ways of observing quantum 
criticality in optical lattices. One is through the analysis of equilibrium atomic density 
profiles, and the other is through dynamics of the sample in the quantum critical regime. 
We start in Sec. [2] with a brief review of the theoretical basis of critical universality in 
a system of ultracold atoms in optical lattices. In Sec. [31 we present progress toward 
observing this critical regime from analyses of in situ atomic density profiles. Also, our 
observations on the mass fiow in optical lattices demonstrate the possibility of observing 
dynamics in the quantum critical regime. 




g = t/u 

Figure 1. Illustration of the phase diagram of the superfluid-to-Mott insulator 
quantum phase transition. The phase diagram is shown as a function of the 
temperature T and coupling constant g = t/U^ where t is the tunneling energy 
and U is the on-site interaction. At zero temperature, the phase transition happens 
at the critical point (T = 0,g = gc)] at finite temperatures, the system shows 
universal behaviors in the V-shaped quantum critical regime. The two graphs of 
atomic distributions in optical lattices illustrate the Mott insulator state (left) and 
the superfluid state (right). Our methods to probe quantum criticality are based on 
studying the universal scaling behaviors by changing the temperature of the system, 
indicated by the vertical arrow, and on studying quantum critical dynamics by ramping 
g across the quantum phase transition, indicated by the horizontal arrow. 
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Figure 2. Illustrations of phase diagrams of the Bose-Hubbard model, (a) Zero- 
temperature phase diagram in the ^/U-t/U plane. Here MI denotes the Mott insulator 
state, SF denotes the superfluid state, and vacuum denotes the state with no atoms 
occupying the lattice sites. Density profiles of trapped gases probe lines in the phase 
diagram along the chemical potential /i direction (with fixed t/U), shown by the dashed 
lines. Points A, B, and C are the specific critical points on the phase boundaries 
referred to in the text. At finite temperatures, these critical points expand to quantum 
critical regions, (b) Finite-temperature phase diagram along the fixed-f/C/ line going 
through the vacuum-to-superfluid transition point A at zero temperature. Here vac. 
denotes the vacuum state and QC denotes the quantum critical region, (c) Finite- 
temperature phase diagram along the fixed-i/f7 line going through the vacuum-to- 
superfluid transition point B and the tip of the occupation number iV = 1 Mott 
insulator lobe (point C) at zero temperature. 



2. Probing quantum criticality in the density profiles and dynamics of 
ultracold atoms 

We explore quantum critical behaviors of ultracold atoms near the superfluid-to-Mott 
insulator transition in a 2D optical lattice. Based on high-resolution in situ absorption 
imaging, we can record density profiles and fluctuations of the sample in equilibrium, 
which allows us to identify critical scaling laws of thermodynamic observables and 
compare to the theoretical predictions. Another direction is to initiate a dynamic 
passage to the quantum critical regime, from which we attempt to extract dynamical 
critical exponents and quantum transport coefficients. In this section, we will discuss 
the "scaling" and the "dynamics" approaches separately. 

2.1. Scaling of equation of state 

Sufficiently near the quantum critical point, thermodynamic quantities, such as the 
atomic density n, obey universal scaling relations [131 El IISl [H] • Consider the equation 
of state n = n{fi,T), where n is the density, /i is the chemical potential, and T is the 
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temperature. Near the superfluid-to-Mott insulator phase boundary, finite temperature 
leads to excitations in the Mott-insulating domain, and the occupation number nd^ 
deviates from the non- negative integer occupation number noc?^ = 0, 1,2, ■■■ in the 
zero-temperature limit. Here d is the optical lattice constant, and the riorf^ = state 
corresponds to the vacuum state. After the non-universal part uq is subtracted from n, 
the remaining part is expected to satisfy the following scaling form in the critical 
regime da M, US] 

T) = T) - riolyU, T) 

where D is the dimensionality, F is a generic function, z is the dynamical exponent, 
u is the correlation length exponent, k-Q is the Boltzmann constant, /ic is the chemical 
potential where the zero-temperature phase transition takes place. The scaling form 
Eq. (JT]) holds when the dimensionality D is below the upper critical dimension D^. = A—z, 
where z will be 1 or 2 depending on what part of the phase boundary the critical point 
resides at. Near the generic vacuum-to-superfluid transition points (see points A and B 
in Fig. El^a)) and the superfluid-to-Mott insulator transition points (other than the tip 
of the Mott insulator lobe), the equations of state are expected to scale according to the 
dilute Bose gas universality class with critical exponents z = 2,u = 1/2; at the tip of 
the Mott insulator lobe (see point C in Fig. |2|^a)), the equation of state scales according 
to the 0(2) rotor universality class with critical exponents z = 1,1^ = 1 [H [131 [13 fT6] . 

We study the critical scaling relations based on in situ absorption imaging of atomic 
density proflles. Applying the local density approximation (LDA) enables deduction of 
the properties of homogeneous systems from the observed properties of non-uniform 
trapped atoms [Ml |T7] . Each density proflle corresponds to a line along the chemical 
potential direction in the phase diagram (for example, see the two red dashed lines in 
Fig. [2]^a)). Sufficiently close to the quantum critical point, when the scaled equation of 
state — ("-"g)'^ — jg plotted versus the chemical potential u for different temperatures, 

all curves should intersect at the same point fi = fic- Once Hc is determined, one can 
plot the scaled equation of state versus the scaled chemical potential StuEAl^^ ^jid all 

curves for different temperatures should collapse into a single curve given by the generic 



function F [T] 

Finally, we point out important issues on observing quantum criticality in optical 
lattice experiments. First of all, since quantum critical scaling laws are valid only 
at sufficiently low temperatures, it remains an interesting question to determine and 
understand the deviation from the expected scaling laws due to finite-temperature effects 
[6l [16] . Recent Quantum Monte Carlo simulations suggest that universal scaling laws 
can be observed at temperatures as high as k^T ~ 6t for atoms in a two-dimensional 
lattice, where the width of the ground band is 8t [16]. Secondly, inhomogeneity of 
the trapping potential can suppress the quantum critical behavior at low temperatures 
when the correlation length is restricted by the length scale of the trapping potential |18j . 
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The finite-temperature and trap effects require the temperature to be small compared 
with the ground band width and large compared with the characteristic energy scale 
of the trapping potential. Finite size scaling algorithms, which have been applied to 
improve the accuracy of analysis in quantum Monte Carlo simulations [19], have also 
been proposed for the experiments [Tl] . 

2.2. Quantum critical dynamics 

We discuss some prospects of studying quantum critical dynamics using cold atoms in 
optical lattices. Prominent dynamic phenomena include quantum critical transport of 
mass and entropy, and dynamics of defect generation across the quantum critical point 
as described by the Kibble-Zurek mechanism [201 |2T] . 
A. Quantum critical transport 

Mass and heat transport across the quantum critical regime provide important 
tests for quantum critical theory p. Sufficiently close to the critical point, one expects 
that transport coefficients obey universal scaling relations independent of microscopic 
physics [H [22]. In two dimensions, in particular, we expect that the static mass 
transport exhibits a universal behavior, in analogous to the prediction on the electrical 
conductivity [231 [21], ^-iid the static mass conductivity at the critical point is given by 
m , , 

. = 1-*.. (2) 

which only depends on the fundamental constants mjli and a dimensionless, universal 
number $cr determined from the universality class of the underlying phase transition. 
Here Ti is the reduced Planck constant, and m is the atomic mass. Analytic predictions 
on the transport coefficients in the quantum critical regime were recently reported on 
the basis of the anti-de Sitter /conformal field theory duality [211 125]. Measurements 
of transport coefficients in general can be of fundamental interest in quantum field 
theory [21]; the relation between mass and thermal conductivities is in close analogy 
to the Wiedemann-Franz relation between charge and thermal transport coefficients in 
electronic systems, which is shown to break down near the quantum critical point in a 
recent experiment [26] . 

Mass and heat transport are induced by generalized forces such as chemical 
potential gradient and temperature gradient. A natural approach to study dynamics of 
atoms in optical lattices is to first create non-equilibrium density distributions in the 
sample and then measure the subsequent evolution of density profiles. 

Non-equilibrium density distributions can be induced in various ways. For example, 
one can create a controlled perturbation in the local chemical potential and induce 
transport by dynamically changing the envelope trapping potential in an equilibrated 
system or changing the on-site interaction IJ near a Feshbach resonance [7]. On the 
other hand, applying lattice ramps slow compared to local microscopic time scales can 
still violate global adiabaticity and induce macroscopic mass and heat flow [2Z] . This is 
aggravated by the pronounced difference in the equilibrium density and entropy profiles 
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between superfluid and Mott insulator phases, as shown in Fig. [3J In a non-equihbrated 
system, we expect quantum critical dynamics to take place near integer site occupation 
numbers. 

While measuring the evolution of the density profile is straightforward using our 
in situ imaging technique [T2l [27], heat or entropy measurement in the quantum 
critical regime remains a challenging task. Nevertheless, the entropy profile is readily 
measurable deeply in the Mott-insulating regime by counting occupancy statistics using 
single-site resolved florescence imaging in combination with on-site number filtering 
[9], [11], or can be extracted from counting average site occupancies before and after on- 
site number filtering processes [27]. Since the local equilibration time scale (on the order 
of h/U jlT]) is sufficiently decoupled from the global dynamics |28], a locally isentropic 
projection from the quantum critical regime deeply into the MI regime can be achieved 
and the local entropy profile measured. 

From the density and entropy profile measurements, we can determine their current 
densities through the application of a generic continuity equation |^ + V- Jp = Tp. 
Here p{x, r) represents experimentally measured mass or entropy density, Jp is the 
corresponding current density, and T p is a source term which characterizes, for example, 
particle loss (Fp < 0) or entropy generation {Tp > 0). 

Mass conductivity a and thermal conductivity n can be determined by relating the 
mass and entropy current densities, Jn and J^, as functions of position x and time r, to 
the generalized forces: the local chemical potential gradient V/i, the potential energy 
gradient W , and the temperature gradient VT. They obey the following transport 
equations [29] 

J„(f , r) = - aV[/i(f, r) + V{x, r)] - ^VT(x, r) (3) 

Kb 

MS, r)= - ^V[/i(x , r) + V{x, r)] - ^VT{x, r). (4) 

Here L„g and Lgn are phenomenological coefficients similar to the Seebeck and Peltier 
coefficients in the thermoelectric effect and can be related via the Onsager reciprocity 
relation [30] . 

Finally, to obtain precise information of spatially resolved chemical potential 
gradient and temperature gradient, we resort to the equilibrium properties of the sample 
which can be determined from measurements of the equilibrium density and density 
fluctuation. The complementary knowledge of the equation of state n{fi,T) and its 
fluctuation 6n'^{fi,T) in equilibrium can be inverted to obtain ^{71,571^) and T{n,6n^). 
We propose that, in a sample driven out of equilibrium globally but remaining locally 
equilibrated, local density and fluctuation measurements can still be used to extract its 
local chemical potential and temperature. This assumption can be further examined by 
comparing local compressibility to density fluctuation and extracting local temperature 
through the application of the fluctuation-dissipation theorem [31j . 
B. The Kibble-Zurek mechanism (KZM) 

Based on general critical scaling arguments, the KZM predicts the formation 
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Figure 3. Sketch of density and entropy profiles of a trapped, finite-temperature 
gas in tlie superfluid regime (SF, blue curves) and the Mott-insulating regime (MI, 
red curves), with the same particle number and total entropy. The gray shaded area 
marks an extended region near unit site occupation number {nd? = 1) where quantum 
critical transport can take place when global adiabaticity breaks down during the 
lattice loading process, d = 532 nm is the lattice constant. 

of topological defects after a system dynamically crosses through a second-order 
thermodynamic [32l |33] or quantum phase transition [20| |2T] . For optical lattice 
experiments, the KZM applies when the system is quenched from a gapped Mott 
insulator state to a gapless superfluid phase, and predicts that the density of defects 
scales with the ramp rate of the coupling strength g |21j . 

The scaling behavior can reveal critical exponents of the underlying quantum phase 
transition [211 EH [35]. When the coupling strength g is adiabatically ramped close to 
gc, the many-body gap A scales as A oc |(7 — gd^'^ = |A(r)|^^, where A(r) = g^r) — gc 
characterizes the time dependence of the ramp. The adiabaticity criterion breaks down 
at a time t = Ta when the gap A becomes small enough and the ramp rate violates 
< A, yielding excitations with a characteristic energy scale A^ oc \\'{Ta)\'^'^^^^'^~^^^ 
or length scale oc The density of defects Uex should therefore scale universally 

as [SHIM] 

ne. cx cx 7^/^^-^+^^ (5) 

where 7^ = |A'(ra)| is the magnitude of the ramp rate at which adiabaticity fails. 

In a two-dimensional superfluid, topological defects are vortices. Observing vortices 
in an optical lattice using in situ imaging is challenging, largely due to the smallness of a 
typical vortex core size (< 1 /im) compared to the imaging resolution (> 5 /im) available 
in most experiments. While the latter can be technically improved, increasing the vortex 
core size by reducing the atomic interaction can also be achieved experimentally either 
through tuning a magnetic Feshbach resonance [7] or releasing atoms for a short time- 
of-flight time [361 EZ]. 

Further extensions of the KZM consider finite-temperature and finite-size effects 
[35] . In general, the scaling of excitations also depends on the pathway of quenching 
[2TI [35I [38] . and the system can enter the Landau-Zener regime in nearly defect-free 
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processes [20]. Detailed experiments could reveal the wealth in the dynamics of quantum 
critical phenomena as well as the intriguing connection between quantum mechanics and 
thermodynamics in genuine quantum systems [391 HQ] . 

3. Experimental progress 

3.1. Scaling of density profiles 

We study quantum critical scaling relations based on in situ imaging of atomic density 
profiles. Our experimental procedure follows three steps. First, we prepare a quantum 
gas of 10^ cesium atoms in a two-dimensional (2D) trap with an aspect ratio of 200 
to 1 by loading a three-dimensional Bose-Einstein condensate into a single site of a 
vertical lattice with a lattice constant of 4 /im. Secondly, we transfer the 2D quantum 
gas into a 2D optical lattice (with lattice constant d = 532 nm) which induces the 
superfluid-to-Mott insulator transition in sufficiently deep lattices. Finally, we perform 
high-resolution absorption imaging to record the atomic density profiles. Details on the 
apparatus and experimental sequence can be found in earlier reports [12], [271 E] • 

Previously, we observed scale invariance and universality in two-dimensional Bose 
gases with different temperatures and interaction strengths [H]. This observation 
motivates us to verify the universal scaling relations near the quantum phase transition 
for atoms in optical lattices. We work with scattering length a = 300aB, lattice depth 
Vhat = S.S-Er where = h x l,326Hz is the recoil energy, and h is Planck's constant. 
At this lattice depth, the energy gap between the ground and first excited bands is much 
larger than the temperature of the system, and the system is governed by ground-band 
physics; we can reach a temperature as low as T = At/k-Q, where t = x 2.7 nK is the 
tunneling parameter. 

Near the vacuum-to-superfiuid transition point A in Fig. [21 the non-universal part 
of the density, no, introduced in Eq. ([T]), equals zero, and the density profile n is expected 
to obey the scaling relation of the dilute Bose gas universality class [T3l [T5] : 



where / is a universal function, the exponents are z = 2, z/ = 1/2. We expect the 
zero-temperature vacuum-to-superfiuid transition in 2D lattices to happen at a critical 
chemical potential fic = ~4t, which can be understood in the following picture. The 
ground band of a 2D lattice has a width of 8t, with the middle of the band chosen as the 
chemical potential zero point. At zero temperature, when the chemical potential is below 
the bottom of the band (/i < — 4t), atoms are not allowed in the lattice and the ground 
state is a vacuum; as soon as the chemical potential reaches — 4i:, atoms are allowed to 
occupy the lowest energy state and will form a superfiuid. Thus the zero-temperature 
phase boundary is fic = ~4t for the vacuum-to-superfiuid transition. At sufficiently low 
temperatures, we expect that the constant-temperature traces of the scaled equation of 
state, when plotted as functions of the chemical potential /i at different temperatures. 




(6) 
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should have a common crossing point at as is suggested by Eq. 

To verify this scahng relation, we take a series of atomic density profiles at different 
temperatures, shown in Fig. Hl^a). We extract the temperature and peak chemical 
potential by fitting the low-density tails according to the following formula, which is 
originally derived for non-interacting gases |12] and is here modified to include the 
mean- field interaction effect, and to include all orders in the fugacity expansion [33]: 

1 oo 

n{x) = E [^0 Wt)f e'^^l^o-^^'^ff^^'^^-^^"^)!, (7) 
^ 1=1 

where n is the 2D atomic density, d = 0.532/xm is the 2D lattice constant, /o(a;) = 
JItt ^^P (3;cos^) is the zeroth order Bessel function with purely imaginary argument, 
/3 = /iQ is the peak chemical potential, Qcs = Ud"^ |33] is the effective interaction 
strength for a 2D thermal gas in a 2D lattice, and V is the trapping potential. 



(a) (b) 




Radius r (d) Normalized chemical potential 



Figure 4. Density profiles and the scaled equation of state at different temperatures 
(at Vlat — 6.8-Er, t/U = 0.18). (a) Occupation number as a function of radius 
measured at different temperatures: T = 29 nK (black circles), 26 nK (red triangles), 
24 nK (green triangles), 18 nK (blue diamonds), 15 nK (cyan triangles), and 11 nK 
(yellow hexagons, with a lower particle number), (b) Scaled equation of state as a 
function of local chemical potential. When fceT w lOt, the constant-temperature 
traces of the scaled equation of state cross at /i w — lOi; when fceT « 5t, the curves 
cross at /z « — shown by the two arrows. 



By plotting the scaled equation of state as a function of the chemical potential /i, 
shown in Fig. Hl^b), we see that the constant-temperature traces of the scaled equation 
of state at different temperatures cross at different chemical potentials. We study 
this temperature dependence by finding the crossing point for each pair of constant- 
temperature traces of the scaled equation of state, and plot the chemical potential at 
the crossing point as a function of the average temperature of the pair, shown in Fig. 
|5](a). As the temperature decreases, the crossing points move toward higher chemical 
potentials. We take the data with the lowest four temperatures, and plot the scaled 
equation of state as a function of the scaled chemical potential ^7^S^'° , where 
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/^cr,o = is the expected zero-temperature transition point, shown by the red dashed 
hne in Fig. |5] (a). Taking /icr,o = '~4t, we can test the scahng law, Eq. (j6]), by rescahng 
the density and chemical potential. At low temperatures (10 nK < T < 24 nK), the 
rescaled data, shown in Fig. [5] (b), are much closer to each other than they are in Fig. H] 
(b), which is consistent with the critical scaling law given by Eq. ([6]), near the critical 
point yUcr,o = ~4t. 

At low temperatures, our preliminary measurement on the equation of state is 
consistent with the expected critical scaling law given by Eq. (E]). To test the critical 
scaling law more quantitatively, it is important to determine the peak chemical potential 
and temperature of the sample accurately in the fitting. Further analysis is underway 
to reduce the systematic errors in the fitting, and to test scaling laws near points B and 
C in the phase diagram (see Fig. [2]). 



(a) (b) 




Normalized temperature k T/t (/^ - /icr,o)/ {k^T) 



Figure 5. Crossing points and overlapping of the constant-temperature traces of 
the scaled equation of state (at Vcat = 6.8£'r, t/U = 0.18). In (a), the chemical 
potential at the crossing points defined in Fig. |4] (b) is plotted as a function of the 
average temperature (black solid circles). The asymptote in the zero-temperature 
limit, /icr,o = ~4t, and the mean-field predictions based on Eq. ([7|) are shown as the 
red dashed line and the blue dashed line, respectively. In (b), the scaled equation 
of state is plotted against the scaled chemical potential based on /icr,o = — 4t for the 
low-temperature data: 11 nK (yellow hexagons), 15 nK (cyan triangles), 18 nK (blue 
diamonds), and 24 nK (green triangles). Here the ground band in a 2D optical lattice 
potential has a width of 8i =21 nK. 



3.2. Dynamics 

Our recent experiment studied global mass transport and statistical evolution in a 
2D sample across the SF-MI phase boundary [27]. We discovered slow equilibration 
dynamics with time scales more than 100 times longer than the microscopic time scales 
for the on-site interaction and tunneling energy. This suggests that transport can limit 
the global equilibration process inside a sample traversing a quantum critical point. 
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Figure 6. Evolution of the density profile and the atom number current density after 
a short 50 ms lattice ramp from zero depth to a final depth of 10 Eji [U jt = 11). Upper 
figure shows the density profile after holding the sample at a constant final depth for 
hold times r =(a) 10, (b) 150, and (c) 350 ms (black circles). In each figure (a-c), 
the ncar-cquilibrated density profile measured at long hold time r = 500 ms (gray 
triangles) is plotted for comparison, (d) shows the atom number current density at 
hold time r = 10 (black squares), 150 (red circles), and 350 (blue triangles) ms, derived 
from the density profiles measured near hold times shown in (a-c) using Eq. ([3]). 



In Fig. [6] (a-c), we plot the evolution of density profiles of a 2D gas containing 
N = 2 X 10^ atoms after a short 50 ms ramp from zero to a final lattice depth of 10 Efj. 
At this lattice depth, U/t = 11 is below the critical point U/t = 17 for the Mott insulator 
state with unit occupation number [19]. We record density profiles after holding the 
sample at the final lattice depth for various hold times r. With an equilibration time 
scale around 180 ms, the cloud gently expands and the peak density slowly decreases due 
to the increase of repulsive atomic interaction during the lattice ramp. This equihbration 
time scale can depend on the sample size and the local properties of the coexisting phases 
in an inhomogeneous sample. 

We further extract the evolution of local mass current density, leading to detailed 
local transport properties beyond a single equihbration time scale. We compute the 
mass current density J„ by comparing density profiles taken at adjacent hold times 
(Ar = 10 ~ 50 ms) and applying the continuity equation, + V- J„(r, r) = 0, 

to evaluate Jn{r,T). Here, we assume no atom loss in the analyses for short hold 
times r < 500 ms. Assuming that mass flow only occurs in the radial direction 
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(f) due to azimuthal symmetry of the sample, we write the mass current density as 



where A^(r, r) = n{r', r)27rr'dr' is the number of atoms located inside a circle of radius 
r at hold time r. Positive j means a current flowing toward larger radius r, and vice 
versa. 

In Fig. E^d), we show j(r, r) computed from density profiles measured near hold 
times shown in Fig. |6] (a-c). We observe overall positive mass flow, which is consistent 
with the picture of an expanding sample inside the optical lattice. The mass current 
density varies across the sample. Shortly after the lattice ramp at r = 10 ms, mass 
transport is most apparent inside a radius r = 40d, where the occupation number 
nd'^ > 1 and the atoms respond to the increase of on-site repulsion. The current density 
peaks around an annular area 20d < r < 30d when the occupation number is in the 
range 2 < nd'^ < 3; outside this annular area, the current density is suppressed when 
the occupation number is in the range nd"^ > 3 or nd"^ < 2. At a larger hold time 
r = 150 ms, similar transport continues to take place but with smaller amplitude. At 
a long hold time r = 350 ms when the sample is closer to equilibration, the current 
density j becomes smaller than our measurement noise. 

In this section, we have shown that spatially resolved mass current density 
is readily measurable using our in situ imaging technique. We expect that local 
transport coefficients can be extracted using Eq. ([3]), from further measurements of 
local temperature gradients and chemical potential gradients. Our interest lies in mass 
transport in the quantum critical regimes near integer occupation numbers, where the 
static mass conductivity is predicted to be universal (Eq. ([2])). Measurements of local 
entropy density are under future investigations, with details outlined in section 12.21 
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